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ABSTRACT 

We present Gran-Telescopio-Canarias/OSIRIS optical spectra of 4 of the most compact and massive 
early-type galaxies in the Groth Strip Survey at redshift z ~ 1, with effective radii R c = 0.5 — 2.4 
kpc and photometric stellar masses M* = 1.2 — 4 x 10 11 M Q . We find these galaxies have velocity 
dispersions a ~ 156 — 236 kms~ . The spectra are well fitted by single stellar population models 
with approximately 1 Gyr of age and solar metallicity. We find that: i) the dynamical masses of these 
galaxies are systematically smaller by a factor of —6 than the published stellar masses using BRIJK 
photometry; ii) when estimating stellar masses as 0.7xMd yn , a combination of passive luminosity 
fading with mass/size growth due to minor mergers can plausibly evolve our objects to match the 
properties of the local population of early-type galaxies. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: kinematics and dynamics - 
galaxies: structure 



1. INTRODUCTION 

The evolution of the structural properties of massive 
early-type galaxies (ETGs) is currently a topic of de- 
bate with many open questions. In recent years, stud- 
ies have shown that galaxies with M± > 10 1:L Mp are 
smaller by a f actor of four at z > 1. 5 (iDaddi e t al. 2005; 
Truiillo et aT1l2~0 06; Lo nghetti et~aLll2007HBuitrago et all 
20081: Damian ov et alll2009T) and a factor of two at z ~ 1 
(|Truullo et all 120071 1201 ID than the nearby population 
of the same mass. Fo r the most compact and massive 
galaxies in particular, ITruiillo et all (|2009f ) found that 
only -0.03% of SDSS galaxies with M,>8x 10 10 M Q at 
z < 0.2 are ETGs with R c < 1.5 kpc. This contrasts 
with the much higher fraction —15% of such galaxies 
found at z~l by Trujillo et al. (2007, hereafter T07). 
Most interestingly, the lu minosity-weighted ag es of the 
local compact ETGs from ITruiillo et all (|2009ft are very 
low (—2 Gyr), which argues against this sample being 
the counterpart of a passively evolved population of high 
redshift galaxies. Thus, this suggests that most high red- 
shift compact and massive galaxies must have undergone 
a systematic growth in size. 

Among the different mechanisms proposed to explain this 
expansion, the one that reproduces best the observed 
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evolu tion of the average mass-size relat ion (ITruiillo et all 
20111) is the minor merger scenario (jNaab et all 120091 : 
Hopkins et all I2010T ). With this mechanism, galaxies 



grow inside-out by accretion of smaller satellites that 
build up an extende d envelope (|van Dokkum et aT1l2010l ; 
iHopkins et all 2009). This model predicts a substantial 
growth in mass and effective radius, with a mild decrease 
in stellar velocity dispersion a. Thus, measurements of 
a at different redshifts are a direct way to constrain this 
evolutionary scenario. In addition, they allow for an in- 
dependent mass estimate through the virial theorem. 
The minor merger scenario is consistent with the com- 
mon findings of a < 250 kms~ for general ETG sam- 
ples at z = 1-2 (iGebhardt et all 120031: iTreu et all 
2005 1 : Ivan der Wei et all [20051 ; Idi Serego Alighieri et all 
20051 iCenarro fc Truiillc l2009t iCappellari et all I200S 
Fernande z Lorenzo et all 



2011) 



Conversely, it is greatly 
challenged by the dramatic evolution required for the 
super-dense galaxies (R c < 1 kpc, M+ > 10 n M m ), which 
are expected to have much larger a. Ivan Dokkum et all 
(2009) measured a ~ 500±100 kms~ from a super-dense 
ETG at z = 2.2, yielding a dynamical mass that is in 
agreement with the photometric mass estimate. How- 
ever, more a measurements of these extreme galaxies 
are needed to support or weaken the need of a stronger 
evolutionary scenario. In this paper we present cr's of 
four massive and compact ETGs at z — 1 , two of which 
are super-dense. We compare dynamical/stellar masses 
and estimate the evolutionary paths of these objects. 
Throughout this work we adopt a standard cosmology 
of Q m = 0.3, n A = 0.7 and H = 70 kms^Mpc- 1 . 

2. THE DATA 

2.1. Sample Selection 

Our targets were selected from the sample by T07, 
which contains photometrically derived parameters for 
831 galaxies with M* > 10 n M© up to z ~ 2 in the EGS 
field from the POWIR/DEEP-2 survey (jBundv et all 
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Fig. 1. — Spectra of our sample galaxies. The black solid lines 
represent the MOVEL fitting results. Light grey regions are [Oil] 
emission at 3727 A and telluric absorption at 7600 A, which were 
excluded from the fits (see text in §2.2). 



2006). Those aut hors measured eff ective radii with 
the GALFIT code (|Peng et al.l I2002T ) and derived to- 
tal st ellar masses by fitting; spec tral energy distributions 
from IBruzual fc Charlotl pJOl with Chabrier IMF on 
BRIJK photometry. We selected four massive and com- 
pact ETGs at z ~ 1 with Sersic index n > 4. Three of 
them have extreme densities in the sense that they are 
the most compact at a given mass. The catalog parame- 
ters of our selection are displayed in the first five columns 
of Table 1. 

2.2. Observations and reduction 

We observed 2 galaxies at a time using optical long- 
slit spectros copy with the OSIRIS camera-spectrograph 
IjCepal [i~99c? ) at the 10.4 m Gran Telescopio Canarias. 
We used a 0.8" slit and R1000 grism in the spectral 
range of 5000 - 10000 A under a seeing of -0.7". The 
total integration time per target was 12800 s. The in- 
strumental resolution was determined by Gaussian fit- 
ting of the strong sky emission lines yielding cri nst = 120 
kms~ . However, in one of the observing nights the 
spectrograph suffered a defocus that decreased the res- 
olution, affecting half of the exposures of 12019899 and 
12024453. For these frames we found <7i ns t = 210 kms -1 , 
and we reduced this data subset independently. The 
standard star L970-30 was also observed to flux calibrate 
the spectra. The data reduction was performed with the 
IRAF longslit package and involved the standard steps 
of bias level correction, spectral flat-fielding, cosmic ray 
removal, sky subtraction, frame co-adding and aperture 
extraction. The telluric absorption band around 7600 A 
was present in all of our spectra, affecting the Ca H+K 
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Fig. 2. — a) Stellar mass-size relation at z ~ 1. Small circles 
represent the entire sample from T07, black squares are our se- 
lection among th em. The straight line is the local relation from 
IShen et al.l 120031) . Black diamonds are our galaxies with stellar 
masses derived as 0.7M dyn . Gray diamonds are their position 
when evolved to z = via minor merger growth, b) B-band Fun- 
damental Plane of ET Gs. The solid line is the local relation from 
Uargensen et all {1996). Our sample is shown as black diamonds. 
Gray diamonds are their position when evolved to z = via mi- 
nor merger growth an d passive luminosity fading. Samples from 
IGebhardt et al.l 112003 1 and Ivan deTWe l et al. (2005) in the range 
z = 0.8 — 1.1 are shown for comparison. We transformed our I 
to B magnitudes appling fc-corrections based on V — I color as in 
Gcbhardt et al. (2003). In addition, we app lied a negative offset o f 
0.1 mag on all surface brightness data from Gcbhardt et al. (2003) 
to transform their Vega to our AB magnitude system. 



lines in targets 12019899 and 12024321. To correct for 
telluric absorption in these two galaxies, we used a F-type 
star observed in the same slit as 12019899. This star's 
spectrum has a signal-to-noise ratio (S/N) = 280 and a 
smooth continuum with no spectral features within the 
telluric band. After its continuum removal and normal- 
ization we obtained the atmospheric transmission profile 
in the spectral window and divided the galaxy spectra 
with affected Ca H+K lines by this transmission pro- 
file. We did not apply the telluric band correction on 
the other targets because the main spectral lines used to 
measure a were not affected in them. Note that this cor- 
rection increases the noise in the spectral window where 
it is applied. Thus, in these cases the correction did 
not contribute towards more precise a measurements and 
for this reason we simply excluded the affected range in 
wavelength from the analysis. The final rest-frame spec- 
tra have (S/N) = 22 — 29 per A and are shown in Fig. 
[TJ Note the strong absorption in the Ca H+K lines and 
the Balmer series from H 7 to H12. [OH1A3727 A doublet 
emission was found in two objects. In 12024790, this 
emission shows a double peak with a separation of 4.1 A 
as compared to 2.7 A of the doublet. 

2.3. Velocity Dispersions 



TABLE 1 
Physical parameters of the observed targets 



Redshift J(ABmag) -R c (kpc) logM*(M ) SB e (mag) logM«(M ) cr(km s" 1 ) logM dyn (M ) Age (Gyr) [Fc/H] 



ID 



12024790 0.9656 

12024321 0.9159 

12019899 0.9325 

12024453 0.9056 



21.89 


0.514 


11.07 


15.82 


10.65±0.03 


156±10 


10.16±0.08 


0.8 


0.0 


21.34 


2.462 


11.10 


18.95 


10.67±0.04 


166±12 


10.90±0.08 


1.0 


0.0 


21.71 


0.469 


11.34 


15.63 


10.80±0.06 


236±17 


10.48±0.08 


1.6 


0.2 


20.91 


1.601 


11.60 


17.60 


10.98±0.05 


186±10 


10.80±0.07 


1.3 


0.2 



Note. — The first five columns are properties from T07, where the effective radii and stellar masses have uncertainties of 0.06 dcx and 0.2 dcx, 
respectively. M^ refers to the stellar mass derived with photometry from the RAINBOW database. SB e is the rest-frame i5-band average surface 
brightness within the effective radius. Age and mctallicity arc luminosity-weighted by the linear combination of SSP models. 



To measure cr's we used the M0VEL code within the 
REDmE distribution (Cardicl 1999). This program com- 
putes an optimal template by creating a linear combi- 
nation of model stellar template spectra that yields the 
best fit to the galaxy spectrum. The linear combination 
coefficients of the templates weigh their contribution to 
the total stellar population and provide an estimate of 
the age and metallicity distribution of the galaxy. Our 
single stellar population (SSP) templates were obtained 
from the MILES libr ary^ (|Sanchez-Blazquez et aLll2006t 
Vazdckis et al. 2010) assuming Chabrier IMF and broad- 
ened to match the resolution of our galaxy spectra. 
These templates were selected with metallicities [Fe/H] = 
{0.0, 0.2} and a range of ages 0.5 — 5 Gyr. We ruled 
out lower metallicities s ince they are very rare given the 
masses of our galaxies (jZahid et al.ll201ll ). 
We fitted each spectrum using a rest-frame range of ap- 
proximately 3700 — 4400 A since the noise generated by 
the subtraction of sky OH lines washed out most fea- 
tures in the redder part of the observed wavelength range. 
Where applicable, we excluded the [Oil] emission and 
uncorrected telluric absorption from the M0VEL fits. In 
order to account for the dominant photon noise in the 
final <j results, we performed Monte Carlo computations 
with 200 bootstrapped spectra. The output produced a 
distribution of values, whose mean and standard devia- 
tion give the value and error of the a. For the objects 
12019899 and 12024453 we calculated a as the weighted 
mean of the values derived from the two different in- 
strumental resolutions, which favors the higher resolution 
value due to its smaller error. T he fits to the spe c tra ar e 
displayed in Fig. Q] Following |j0rgensen et al.l (|1995l ). 
we apply an aperture correction to a diameter of 1.19 
h^ 1 kpc, which increases our values around 6%. The fi- 
nal cr's are given in Table Q] We note that there exists 
publi cly available DE EP2 spectra for all of our galax- 
ies (jDavis et al.ll2007l ). The target 12024453 has already 
been part of the photome t ric an d spectroscopic study by 
Fernandez Lorenzo et al.l (|2011f ). who find a = 173 ± 16 
kms~ from a spectrum of (S/N) = 11 per A. This is con- 
sistent with our measurement. The rest of the DEEP 2 
spectra of galaxies in common with our sample have 
(S/N)< 9 per A and cannot be used for reliable a mea- 
surements. 

2.4. Stellar Populations 

The age and metallicity of each galaxy were determined 
from the values of the SSP model with the largest coeffi- 
cient in the linear combination that best fitted the galaxy 



spectrum (see values in Table 1). These values show rel- 
atively low ages in the range of 0.8— 1.6 Gyr, with uncer- 
tainties estimated to be around 0.4 Gyr. In addition, we 
calculated photometric stell ar masses using the data from 
the RAINBOW databass£3 (jBarro et al.ll201lD . which in- 
cludes photometry for our galaxies in ^30 bands cover- 
ing a range from 150 nm up to 70 /zm. We fit these 
data t o stellar population models bv lBruzual fc Ch ariot 
(2003) with Chabrier IMF , assuming an extinction law 
from lCalzetti et ail (|2000f h We constrain the metallicity 
in the same way than for the M0VEL templates. The re- 
sults show solar metallicity for all galaxies and ages in 
the range of 0.7-1.7 Gyr. Such outcome is very similar 
to the results from M0VEL. 

The galaxies 12024790 and 12024321 have [OII]A3727 A 
emission, implying ongoing or recent star formation. We 
combine the [Oil] fluxes from our spectra with the rest- 
frame ^-corrected 24 /mi fluxes from RAINBOW to de- 
rive star formation rate s (SFR) using the conversion from 
IKennicutt et all (|2009D based on Kroupa IMF. The total 
SFR[oii]+ir for the [Oil] emission galaxies is 9.7 and 22.6 
M0yr -1 , respectively. 

3. DISCUSSION 

Dynamical masses are calculated under the assumption 
of homology as Md yn = /3R e a 2 /G wit h (3 = 5 (see Ta- 
ble 1 ), following studies of local ETGs (Cappcllar i et al.l 
l2006f ). As noted in § 2.2, the galaxy 12024790 shows a 
double peak in the [Oil] emission. If such feature is due 
to rotation of a gaseous disk, it yields a projected veloc- 
ity of V r sin(i) ~ 165 kms~ , which is comparable to its 
measured a = 156 kms~ . However, it is unlikely that 
this caused an important bias in the dynamical mass es- 
timates. At z = it has been fou nd that /3 = 5 is robust 
again st rotation up to V r /a « 1 (|Cappellari et al.ll2006l 
l2007f ) , and there is no clear evidence that this would b e 
different at z = 1 (|van der Wei fc van der Marell [20081) . 
Therefore, we assume that any deviations from homology 
due to possible rotational support are within the errors 
of our dynamical masses. 

Our galaxies have values that are systematically smaller 
than the originally published stellar masses by an average 
factor of ~6. This is much larger than the combined error 
of photometric stellar mass and radius (typically a factor 
of ^3). Note that for the galaxy 12024321 these masses 
are consistent within uncertainties. The stellar popula- 
tion models we used with the RAINBOW data yield stel- 
lar masses in the range of logM* = 10.65 — 10.98, on av- 
erage a factor of ~1.8 larger than the dynamical masses. 



9 miles. iac.c 



rainbowx.fis.ucm.es 



The masses from RAINBOW are derived with the same 
IMF and stellar population models as in T07 but with a 
broader spectral range (UV-FIR) as described in § 2.4. 
In order to investigate the large differences between both 
sets of masses, we also derived masses using only BRIJK 
bands as in T07. We found no significant difference be- 
tween our UV-FIR and BRIJK masses, implying that 
the discrepancy with T07 cannot be simply explained by 
the use of a different spectral range. We also checked 
that our mass determinations are weakly sensitive to the 
choice of extinction mode l and small variations in metal- 
licity. iBarro et al.l (|2011( 1 derive masses for hundreds of 
galaxies in T07's sample using the full RAINBOW pho- 
tometry. They find that, on average, their masses are 
consistent with those from T07. Therefore, the cause 
of the inconsistency with our photometric masses is un- 
clear. Given the large uncertainties in the systematics 
of photometric masses, we adopt dynamical masses as 
the corr ect total mass estima tes. Following local lensing 
studies (jGavazzi et al.l I2007T) . we apply a fiducial dark 
matter fraction /dm = 0.3 and recalculate the stellar 
masses as 0.7 x Mdyn- 

Comparison studies between stellar and dynamical mass 
of less dense ETG s than those in our sam ple have been 
condu cted at low dCappellari et al.l [2006) and high red- 
shifts (jvan der Wei et alJl2006H with the result of dynam- 
ical masses being on average equal or larger than stel- 
lar masses, as expected when considering a dar k mat- 
ter contribution. When Ivan der Wei et all (|2006f l calcu- 
late stellar masses using similar multi-band photometry 
and stellar population parameters as in T07, they find 
M*/M dy n = 0.66 (Salpeter IMF). Therefore, our results 
point towards a systematic overestimation of the stellar 
masses used to select the galaxies of our sample. We 
cannot rule out a deviation from homology that requires 
f3 ^ 5. However, the theoretical expect ation for objects 
with such high Sersic i ndex is /3 < 5 (|Cappellari et al.l 
l2006t iBertin et al.ll2002ll , which would yield even smaller 
dynamical masses. 

Fig. [2^ shows the position of our objects in the stel- 
lar mass-size relation with the old/new stellar mass es- 
timates. Fig. [2}d shows our objects in the Fundamen- 
tal Plane (FP). At a given R c , they have systemati- 
cally larger masses and smaller combination of surface 
brightness (SB e ) and a than the respective local rela- 
tions. Moreover, our objects follow a trend in which 
they lie further from both local relations at smaller radii. 
This result has also been fo und in z ~ 0.9 clusters by 
Uorgensen et al.l (|2006l 120071) . Compared to the FP de- 
rived by these authors, at a given R c our field sam- 
ple shows a similar slope and a relative offset of ~ 0.3 
mag towards brighter SB e , which is in agreement with 
field galaxies having generally younger stellar popula- 
tions than clusters of the same redshift. 
In order to understand how our galaxies relate to low 



redshift ETGs, we investigate the simple evolutionary 
scenario based on passive luminosity fading of the stel- 
lar populations and growth via minor mergers that has 
been pro posed for the general population of ET Gs at 
z ~ 1 (|Gebhardt et alJ 12003b iNaab et all 12009( 1. To 
derive the amount of fa ding for our sample, we use 
iBruzual fc Charlotl (|2003t l models with Chabrier IMF 
and the redshift, ages and metallicities from Table [TJ 
Passive evolution from z ~ 1 to z = results in an aver- 
age dimming of ^2 mag in SB e . We assume that the pas- 
sive luminosity evolution does not affect R c . Minor merg- 
ers are very efficient growing the size at a given mass in- 
crement ,sljghjl^_decreasing (j. According to the simula- 
tions bv lNaab et all (|2009l l. an ETG with M* = IO n Af 
and R c = 1.5 kpc would grow a factor of 1.5 in mass via 
minor mergers from z ~ 1 to z = 0. We adopt this fac- 
tor of mass growth for our galaxies. Following these au- 
thors, we model each merger as cri/a? — (1 + ?/e)/(l + rj) 

and R c j: / R ci = (1 + f/) 2 /(l + r?e), where the subindexes 
(i)/(f) denote initial/final values of the host galaxy and 
e = er^/of and the mass ratio r\ — M a /Mi take into 
account the accreted system (M a ,a a ). We assume a his- 
tory of four mergers with constant rj = 0.1 and e = 0.2 
for all galaxies (note that r\ and the total mass growth 
set the total number of mergers) . T his choice o f values 
is consistent with the simulations bv INaab et al.l |2009). 
The total change for R B and a become factors of 2 and 
0.8, respectively. These results are robust against vari- 
ations in the values of r\ and e at the given total mass 
growth. We assume that the mass-to-light ratio does 
not change significantly after each merger. In addition, 
we assume that as t he size increases, SB P decrea ses by 
ASB e oc 2.94Alogi? e (|Hamabe fc Kormendvlll987l l. 
Fig. [3 shows the combined evolution of minor mergers 
and luminosity fading for our galaxies. Most of them 
reach positions that are consistent with the local FP and 
mass-size relations, considering the estimated uncertain- 
ties in our measurements. Therefore, this simple evolu- 
tionary scenario can plausibly describe the evolution in 
our sample galaxies to become normal ETGs at z = 0. 
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